Mira's wind explored in scattering infrared CO lines 
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ABSTRACT 

We have observed the intermediate regions of the circumstellar envelope of Mira (o Ceti) in 
photospheric hght scattered by three vibration-rotation transitions of the fundamental band of 
CO, from low-excited rotational levels of the ground vibrational state, at an angular distance 
of /3 ^ 2" — 7" away from the star. The data were obtained with the Phoenix spectrometer 
mounted on the 4 m Mayall telescope at Kitt Peak. The spatial resolution is approximately 
0.5" and seeing limited. Our observations provide absolute fluxes, leading to an independent 
new estimate of the mass- loss rate of approximately 3 x 10~^ Mq yr~^, as derived from a simple 
analytic wind model. We find that the scattered intensity from the wind of Mira for 2" < (3 < 7" 

decreases as which suggests a time constant mass- loss rate, when averaged over 100 years, 
over the past 1200 years. 

Subject headings: stars: AGB and post-AGB — circumstellar matter — stars: individual (o Ceti) — 
stars: late-type — stars: mass loss — infrared: stars 



1. INTRODUCTION 

Ever since the circumstellar lines of the a Her 
M supergiant were identified in absorption against 
the photospheric continuum of the G III secondary 
(Deutsch 1956), there have been numerous at- 
tempts to map circumstellar regions. Three spec- 
troscopic tools have been used extensively. The 
most extensively employed are microwave emis- 
sion lines. The IR continuum from circumstellar 
dust, especially in the observations by the IRAS 
and ISO satellites, has been used. Finally, ap- 
proximately a dozen circumstellar shells have been 
imaged in photospheric light scattered by atomic 
resonance lines. In this paper we shall extend the 
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latter technique to resonance scattering in molecu- 
lar lines using the mid-infrared vibration-rotation 
lines of the ubiquitous CO molecule. 

In his review of circumstellar envelopes and 
AGB stars, Habing (1996) noted both the useful- 
ness and the uniqueness of measuring scattered 
resonance lines, for example of Nai and Ki, but 
drew attention to the relatively few cases where 
this method has been applied. This is certainly 
because the observations and their interpretation 
are quite complex. There are uncertainties in re- 
lating these observations to the mass and struc- 
ture of the expanding gas since atoms are depleted 
by dust grains and the details of ionization in the 
circumstellar shell are not known. An advantage 
of microwave measurements is that the molecular 
species sampled may be much less affected by these 
uncertainties. This advantage is also present in the 
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study of photospheric light scattered by circuni- 
stellar CO molecules, as presented in this paper. 
The spectral energy distributions (SEDs) of AGB 
stars have a maximum in the infrared which also 
suggests infrared CO observations. The infrared 
vibration-rotation lines of CO, as an alternative 
to rotational CO lines at millimeter wavelengths, 
also allow the study of regions close to the star and 
admit higher spatial resolution in single telescope 
studies. 

Pioneering investigations of atomic resonance 
line scattering were made of the Kl 7699A circum- 
stellar emission in a Ori by Bcrnat & Lambert 
(1975, 1976). Recent work includes the study of 
circumstellar shells around three N-type carbon 
stars using Ki 7699 A (Gustafsson et al. 1997) and 
of circumstellar shells of M-type mira stars using 
various atomic resonance lines (Plez & Lambert 
1994; Guilain & Mauron 1996). 

Imaging off-star emission of CO at 4.6 /im has 
previously only been done a few times (Sahai & 
Wannier 1985; Ryde et al. 1999). Dyck et al. 
(1983) performed speckle interferometry on a 
number of CO fundamental vibrational-rotational 
lines. Sahai & Wannier studied the intermedi- 
ate regions of the circumstellar envelope of the 
dust-enshrouded and bright IR star CW Leonis 
(IRC+10 216) by using an annular aperture, and 
determined a kinetic temperature of its shell at 
an average radius of 2". This technique provided 
only spatially averaged fluxes but led the authors 
also to the conclusion that the mass-loss rate in 
the inner parts of the circumstellar shell is less 
than that corresponding to the region observed in 
millimeter CO lines. Since the Sahai and Wan- 
nier observations, great advances have been made 
in infrared detectors. It is, however, not trivial 
to apply infrared arrays in high resolution spec- 
troscopy since at wavelengths longer than about 
1.6 yum thermal radiation from room tempera- 
ture spectrometers dominates the stellar signals. 
Cryogenic spectrometers are therefore required, 
such as the Phoenix spectrometer (Hinkle et al. 
1998). One of the justifications for building this 
spectrometer was, in fact, to map circumstellar 
shells as presented in this paper. In Ryde et al. 
(1999) we demonstrated methods of observation 
and analysis of the vibration-rotation CO emission 
lines from the intermediate regions of circumstellar 
winds. Here, we present results of an application 



of these methods to the circumstellar envelope of 
the M4-7IIIe giant o Ceti (Mira). 

Section 2 discusses Mira. Sections 3 and 4 de- 
scribe the observational set-up and the reduction 
procedures, and the absolute flux calibration, re- 
spectively. Section 5 reviews the observational re- 
sults and Section 6 discusses these results in the 
light of an analytic model for circumstellar en- 
velopes. 

2. MIRA 

A study of the circumstellar envelope of o Ceti 
(Mira or HD 14386) is of interest for several rea- 
sons. The star is the prototype of the mira-class 
of long-period variables characterized by large am- 
plitude visual variations. In the case of o Ceti 
the V-band amplitude is more than 6 magnitudes 
and the mean period is 332 days. The miras are 
cool AGB stars with most of the energy emitted 
in the infrared. The brightness of Mira reflects its 
closeness, (128 ± 18) pc (ESA 1997). The stellar 
mass loss produces a circumstellar envelope which, 
at the distance of Mira, has an angular extension 
of 2' on the sky (Loup et al. 1993). Note, how- 
ever that o Ceti departs from many miras in be- 
ing a binary system. The angular distance be- 
tween Mira A and the hot, compact companion 
star Mira B (VZ Ceti) is 0.6" (Karovska et al. 
1993; 1997). The companion could be a white 
dwarf with a mass of about 1 Mq , a luminosity 
of 2 L0, and a temperature of more than 30 000 K, 
which is embedded in an accretion disk, giving rise 
to an abnormal illumination of Mira A (Danchi 
et al. 1994). Also, note that a number of other 
high-resolution, optical and infrared wavelength 
measurements provide evidence that Mira itself 
is elongated (see discussion by Ryde & Schoier 
(2000) and references therein). 

Theoretical model calculations (Bowen 1988; 
Hofner & Dorfi 1997; Hofner et al. 1998; see 
Woitke (1998) for a review) indicate that the mass 
loss of these stars is most probably caused by 
a combination of radial pulsation and radiation 
pressure on dust grains and/or molecules. It may 
seem natural to expect this mass loss to be spher- 
ically symmetric, but stellar rotation may intro- 
duce a latitude dependence (Dorfi & Hofner 1996). 
Also, on the scale of the spatial resolution of 
our observations, corresponding to a time-scale of 
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about 100 years, one could expect the wind to be 
homogeneous (observations of Mira show the cur- 
rent pulsational behavior dates back to at least 
1638). 

However, there are indications from radio ob- 
servations that the wind of Mira departs from this 
simple picture. There is an asymmetry and a 
shoulder clearly visible in the radio-line profiles, 
see for example Planesas et al. (1990a). There 
have been several suggestions in the literature of 
different multi-wind scenarios or other phenom- 
ena introducing dramatic variations in the enve- 
lope structure with distance in order to fit the 
profiles. There exist several combinations of mass- 
loss rate, expansion and turbulent velocities that 
are able to reproduce the Mira circumstellar line 
profiles rather well and as a result there is a lack 
of consensus regarding the actual expansion ve- 
locities of Mira's wind. All these results are based 
on microwave and sub-millimeter observations of 
rotational CO lines. Crosas et al. (1997) exper- 
iment with Vexp ~ 2 and vturb ~ 4kms~-^ for 
their inner wind, while Young (1995) arrives at 
Vexp = 4.8kms~^ from a fit. However, he points 
out that the wings show an expansion velocity of 
lOkms"^. Knapp et al. (1998) model the cir- 
cumstellar envelope (CSE) with a fast outer wind 
with a mass-loss rate resembling single wind miras, 
M = 4.4-10"^ M0 yr~^, and an expansion velocity 
of (6.7±1.0) km s-^ The inner wind is supposed 
to be a resumed wind in analogy with the detached 
shells found in four carbon stars (Olofsson et al. 
1996). This slow wind component has a lower 
mass-loss rate and a much lower expansion veloc- 
ity than what is found normally in miras; M = 
9.4-10-8 Mq yr-i and Vexp = (2.4±0.4) km s'^ 
Planesas et al. (1990a) interpret their observa- 
tions of CO(J =1-0) and C0( J = 2-1) lines as 
partly originating from a dominating, spherically 
symmetric CSE with Vcxp = 3 km s^^ and a partial 
collimation of the stellar wind into an additional 
bipolar lobe. 

Obviously it would be of great interest to 
study whether the peculiarities observed in the 
microwave lines are also reflected in the infrared 
CO vibration- rotation lines, which map the inter- 
mediate regions (approximately 100 — 1000 stellar 
radii) of Mira's circumstellar environment. 

The geocentric, radial velocity of Mira is large 
enough for the stellar CO lines to be shifted out 



of the telluric lines. This is a key factor since the 
blending of the telluric and stellar lines makes the 
analysis of the faint circumstellar emission very 
difficult. 

Multi-wavelength angular diameter measure- 
ments of Mira (Haniff et al. 1995) combined with 
Hipparcos trigonometric parallaxes (Van Leeuwen 
et al. 1997) suggest a mean radius of (464±80) Rq. 
Based on the work of Mahler et al. (1997) on ra- 
dius and luminosity variations of Mira from Wing 
near-IR photometry, one can derive a bolometric 
luminosity Ltot = 8900 L©, a temperature Tgg = 
2400 K, and a radius of 550 Rg = 3.8 • 10" cm for 
the phase of o Ceti on the date our observations 
were made, 1998 October 28-29. These values are 
consistent with the corresponding values of, for 
example, Danchi et al. (1994) and Haniff et al. 
(1995). They arrive at similar temperatures; a 
mean (over phases) effective temperature of about 
Teff ~ 2800K. We note, however, that this value 
is low as compared with the values suggested by 
Perrin et al. (1998) for giants in the interval M4HI- 
M7ni which is characteristic of miras. Assuming 
a mass of « 1.0 M©, we find the logarithmic 
surface gravity, log 5 = log(GM/R^), to be ap- 
proximately — 1.0 [cgs]. 

3. OBSERVATIONAL SET-UP AND RE- 
DUCTIONS 

The observations were carried out in 1998 on 
October 28 and 29 using the 4 meter Mayall tele- 
scope at Kitt Peak with the Phoenix spectrometer, 
mounted at the Cassegrain focus of the equato- 
rially mounted telescope. This cryogenic, single- 
order echelle spectrometer is a long-slit, high spec- 
tral resolution instrument, designed for the 1 — 
5 /tim region (Hinkle et al. 1998) and marks a ma- 
jor achievement. At 4.6 /xm it is now possible to 
observe objects 5 magnitudes fainter than was pos- 
sible with the Fourier Transform Spectrometer on 
the same telescope. 

The detector is an Aladdin 512x1024 element 
InSb array cooled to 35 K. The rest of the spec- 
trometer is cooled to 50 K. For our observations, a 
30" long and 0.4" (2 pixels) wide entrance slit was 
used, resulting in a spectral resolution of around 
R = 60 000. The slit width projection onto the 
sky of 0.4" means that the spectral information 
is an average over, at least, a region of this size. 
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The visual seeing was around 0.8" , but we note 
that the seeing disk in the infrared (4.6 /zm) is 
60% of that in the visual (Lena 1988), so that the 
seeing nearly matched the slit size, which is ap- 
proximately 30% greater than the diffraction limit 
of the telescope. The dispersion of the spectro- 
graph leads to an over-sampled observed spectrum 
and each spectrum covers a very small wavelength 
range; approximately 11 cm~^. 

We observed the circumstellar envelope of 
Ceti in a hashed configuration, resembling a 
perpendicular railway crossing with the star in 
the middle, see Figure 1. The long slit was placed 
in the off-star position at 2" away from the star. 
We were able to detect emission to a position cor- 
responding to a maximum distance from the star 
of 7" . The four positions will overlap at a distance 
of approximately \/2'^ + 2^" away from the star, 
which enables us to make a relative calibration 
of the observations, for the different slit positions 
around the star. 

We chose an echelle setting for observing low 
excitation lines of the vibration-rotation, funda- 
mental R-branch of the electronic ground state 
of The lines selected, 1-0 R(l) (2150.86 

cm-\ i.e. v= 1 ^ and J' = 2 ^ J" = 1), 
1-0 R(2) (2154.60 cm^i), and 1-0 R(3) (2158.30 
cm~^), have minimal interference with telluric 
lines. The potentially important 1-0 R(0) line has 
an interfering telluric water line. Thus, we ob- 
served the region 2150 to 2160 cm~^, the spec- 
tral coverage being limited by the array length. 
The spectrometer works in order 12 at these wave- 
lengths. 

The CO fundamental bands are located in the 
thermal infrared where the sky background radi- 
ation is high and somewhat variable. Phoenix 
is not a sky balanced device, which means that 
the thermal radiation is recorded. At the 4.6 /um 
wavelength of CO, the thermal background ra- 
diation at Kitt Peak would limit exposure times 
at R = 60 000 to about ten minutes. However, 
the 4.6 /um spectrum has a contribution of telluric 
lines, mostly CO and H2O, which are seen in emis- 
sion. To avoid saturation of the telluric lines the 
exposure time was limited to 60 seconds. 

Our observations consist of two types: on-star 
and off-star exposures. On-star exposures are re- 
quired to remove scattered light from the off-star 
exposures. In order to remove the thermal spec- 



trum from the source spectrum, two exposures 
are required, one of the source and one of the 
sky (Joyce 1992). All spectra must also have ar- 
ray pixel sensitivity differences removed using flats 
and darks as described by Joyce (1992). Since 
Phoenix has a long slit, nodding along the slit 
is possible for point sources. Therefore, on-star 
spectra were taken by moving the star along the 
slit by 15" between exposures. The brightness 
of Mira limited these exposures to 30 s. Each 
set of exposures then gave, after differencing, two 
background-subtracted, on-star spectra. Since the 
stellar signal is cxtincted by the optical depth of 
the telluric lines, the background-subtracted spec- 
tra show in absorption the telluric lines which ap- 
pear in emission in the raw images. 

Off-star frames were taken the first night by 
alternatively observing 2" south and making an 
identical background observation 62" south of the 
star, both with an exposure time of 60 s. Ceti 
has a CO shell with a radius of 57", as measured 
by radio observations (Loup et al. 1993). The total 
effective observing time in the south position was 
66 minutes. Corresponding observations were sub- 
sequently made in the north position, 2" N, also 
here with an exposure time per nod of 60 s. The 
total effective observing time in the north position 
was 80 minutes. The second night a more efficient 
nodding algorithm was employed. Off-star frames 
were taken by alternatively observing 2" east and 
2" west, 60 s each, then a background observation 
was taken 60" west off the star. A total of 54 effec- 
tive minutes were observed in both the west and 
east positions. Over the few minute (at maximum) 
interval between the source and sky observations 
the sky background level is almost constant and 
the airmass changes only slightly. While the ther- 
mal radiation and extinction of well- mixed telluric 
gases, such as CO, are scaled with airmass and 
time, water is not, on most nights being in the 
form of clouds visible in the infrared. As a re- 
sult the telluric water vapor lines will not cancel 
cleanly in the subsequent differencing of spectra. 

The on-star spectra were used to remove the 
effects of scattered stellar light in the terrestrial 
atmosphere and in the spectrometer. Thus, these 
were normalized to unit level and then scaled to 
the levels of the off-star spectra. The circumstellar 
CO vibration-rotation emission was obtained by 
subtracting the scaled on-star spectra from the off- 
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star ones. The emission we wish to detect is about 
50 times weaker than the on-star intensity. The 
on-star intensity is typically about twenty to thirty 
times brighter than the sky background. 

The reduction of the data was performed using 
standard routines of the latest version of IRAF 
(Image Reduction and Analysis Facility). The 
wavelength calibration was made by using the tel- 
luric lines in every frame. The accuracy of the 
wavclcngtli calibration is of the same order as the 
resolution. Therefore, the uncertainty in the wave- 
length scale is probably approximately 0.04 cm~^ 
or 6 kms~^. 

From the ISO data archive^ we have retrieved 
the spectrum around and the flux at 4.6 /um of 
Mira, obtained with the Infrared Space Observa- 
tory (Kessler et al. 1996). The reductions were 
made using the most recent pipeline basic reduc- 
tion package OLP (v. 7) and the ISO Spectral 
Analysis Package (ISAP v. 1.5). The pipeline pro- 
cessing of the data, such as the flux calibration, 
is described by de Graauw et al. (2000); the com- 
bined absolute and systematic uncertainties in the 
fluxes are of the order of ±10%. 

4. ABSOLUTE FLUX CALIBRATION 

The absolute calibration of the spectra in the 
4 fj,m. region was made by a comparison of the mea- 
sured on-star flux (in counts) with the absolute 
flux (in physical units) expected in the broad M- 
band. The conversion factors obtained are then 
used for calibrating the off-star spectra. This 
method contains a number of uncertain steps, re- 
lated to the variability of Mira, the uncertainties 
in the calibration of broad-band photometry, the 
scaling from the broad-band to the narrow wave- 
length region observed spectroscopically, the ques- 
tion of how much stellar light was caught within 
the spectrograph slit in the on-star observation, 
and the question whether the sensitivity of the 
detector or the transmission of the Earth's atmo- 
sphere might have changed between the on-star 
and off-star observations. In order to check the 
calibration we therefore compared with ISO fluxes 
and made a final test in a comparison between the 

^Partly based on observations with ISO, an ESA 
project with instruments funded by ESA Member 
States and with the participation of ISAS and NASA, 
http; / /isowww.estec.esa.nl/ 



observed spectrum and a calculated one from a 
model atmosphere, allowing for the stellar radius 
and the distance as measured by the Hipparcos 
satellite. 

The amplitude of the Mira light curve is 
visually large, but decreases rapidly towards 
longer wavelengths. From the NASA Catalog 
of Infrared Observations (Gezari et al. 1984) 
a mean of eleven measurements gives a flux 
in the Johnson M filter of logj^i, = 3.6 ± 0.2 
[Jy]. This M-band mean fitix is easily con- 
verted to {5.6'tl]f) X 10^^ ergs"-^ cm~^ /xm~-^ or, 
since 1 pixel corresponds to 0.26 A, (1.5^q'5) x 
10^^^ crgs^^ cm^^ pixel^^. In view of the fact 
that the fiux is peaked towards the blue end 
of the M-band in such a way that the fiux in 
the spectroscopic band is higher than the mean 
flux in the band (cf. Figure 2), one must cor- 
rect the absolute calibration. From the low- 
resolution ISO spectrum of this region of Mira 
we find a correction factor of 1.23. We ob- 
tained this correction by computing a mean of 
the ISO spectrum over the M-band and compar- 
ing with the ISO flux at 4.64 /im. Our model 
atmosphere (see discussion below) gives a similar 
result. Thus, for October 28 we arrive at a con- 
version factor of (3.3^^:^) x 10"" erg cm ^ per 
detector count and for October 29 we get a factor 
erg cm ^ per detector count. 
This leads to a fiux in the 4.64 /zm band of 4900 Jy 
or 6.9 X 10"^ erg s"^ cm~^/xm~^. This may read- 
ily be compared with the ISO flux at the same 
wavelength, (4700 ± 400) Jy, cf. Section 3. 

In order to check the calibration further we also 
used an OS-MARCS spherical model-atmosphere for 
the stellar parameters Tes = 2400 K, log 5 = — 1.0 
[cgs], and solar metallicity and allowing for the 
stellar radius (R* = 3.8 • 10"'^^ cm) and a Hippar- 
cos distance d = 128 pc, to generate a synthetic 
spectrum of the region of interest. The stellar pa- 
rameters are determined for the phase of o Ceti at 
October 29, 1998, the date the NIR observations 
were made, see Ryde & Schoier (2000). The model 
atmosphere is a part of a new grid of model atmo- 
spheres (Plez et al. 2000) being generated by an 
extensive update of the marcs code and its input 
data (based on Gustafsson et al. 1975). Molecu- 
lar lines of II2O (with line lists of Partridge and 
Schwencke 1997), of CO (lists of Goorvitch 1994), 
and of many other molecules, are taken into ac- 
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count in the spherical radiative transfer for the 
calculation of a synthetic spectrum. 

Figure 3 shows the spectrum generated from 
the model atmosphere and a comparison with the 
observed one. The general features are well re- 
produced and are mainly due to photospheric wa- 
ter vapor and CO. The agreement of the over-all 
flux level is astonishingly good; the model flux is 
around 30% too small. We find this agreement be- 
tween the empirical, calibrated, absolute flux and 
the flux we calculate from a model atmosphere 
very satisfactory, in view of all the uncertainties 
anticipated. 

It is worth noting that synthetic spectra from 
the new generation of spherical MARCS models in 
the near infrared are able to reproduce observed 
spectra fairly well, even for a pulsating star like 
Mira. While the cyclic variations in mira spectra 
are well known and conspicuous at high resolution 
(Hinkle et al. 1982; Hinkle et al. 1984), little is 
know about phase dependent variations of the 4.6 
fim photospheric spectra of mira stars. It is pos- 
sible that depths of formation of the continuum 
and lines at 4.6 /um are in a region of the photo- 
sphere similar to that of the visual spectrum where 
the dynamic behavior does not have large effects 
on absorption line formation. However, it is also 
possible that the phase of observation gave a fortu- 
itous match to the model. In either case the stellar 
atmosphere is extended and the MARCS match to 
the spectrum is impressive. 

Figure 4 shows a blow-up of the photospheric 
CO R(2) absorption line as well as the scattered 
circumstellar R(2) emission. This emission line 
is shifted by 0.42 cm^^ compared to the telluric 
wavelength scale, corresponding to a radial veloc- 
ity of 58.6kms~^. Mira has a vlsr = 47kms~^ 
deduced from radio data (Loup ct al. 1993), which 
equals a heliocentric velocity of 57kms~^. At 
the time of the observations and at Kitt Peak, 
this corresponds to a geocentric radial velocity of 
SQkms""'^, in excellent agreement with the obser- 
vations. Thus, the wavelength of the scattered 
light is centered on the laboratory wavelength cor- 
rected for the stellar radial velocity. The photo- 
spheric absorption lines, however, are also shifted 
due to the velocity of the pulsating photosphere. 
The bisector of the photospheric absorption line 
lies less than 2kms~^ blue- wards of the center 
of the emission line. The weak emission in the 



absorption line in the on-star spectrum is due to 
scattered light from the CSE. This light is included 
in the on-star measurement since the long-slit will 
also cover the CSE in two directions away from the 
star. Also present in the on-star spectrum is ad- 
ditional absorption along the line-of-sight through 
the circumstellar shell. In our measured, off-star, 
scattered CO light we do not correct for these com- 
ponents. However, they are of nearly equal inten- 
sity and should cancel at the level of the uncer- 
tainties of our measurements. 

Note that when discussing relative fluxes, 
such as in Section 5 where the R(l)/R(3) and 
R(2)/R(3) ratios are calculated, the details in 
the flat-fielding may introduce uncertainties. The 
three lines span over the entire range of the detec- 
tor, and therefore variable effects not taken care 
of by the flat-fielding will show up as an error in 
the ratios. 

5. OBSERVATIONAL RESULTS 

Figure 5 shows the off-star spectrum of the 
west position. Superimposed is the correspond- 
ing, scaled, on-star spectrum, which at least ap- 
proximately represents the radiation exciting the 
molecules. The intensity of this spectrum is scaled 
in order to fit the general features in the off- 
star spectrum. The off-star spectrum obviously 
consists partly of stellar light that is scattered 
in the Earth's atmosphere, in the telescope, the 
spectrometer, and/or by dust grains in the cir- 
cumstellar shell, leading to a spectrum resembling 
the on-star one. Nearly all features are identi- 
fied as photospheric CO and H2O fines, cf. Sec- 
tion 4 and Figure 3. The on-star spectrum in- 
cludes tens of CO vibration-rotation lines of vari- 
ous excited vibrational states. The cold off-star 
spectrum, however, includes only the 1-0 R(l), 
R(2), and R(3) vibration- rotation emission lines 
of ^^CO within the observed spectral range, ex- 
cept for the scattered, on-star light. Thus, the 
circumstellar molecules are radiatively excited by 
the stellar light, which is re-emitted as emission. 

Figure 6 shows the resulting CO emission from 
the observations west, cast, north and south of 
o Ceti, integrated over the long slit. This emis- 
sion is recovered from the data by subtracting the 
scaled on-star spectrum from the off-star one. The 
resulting emission lines are the circumstellar R(l), 
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R(2) and R(3) lines. Variations in the telluric lines 
during the time between the object frame expo- 
sure and the background exposure and/or during 
the time between the on-star and ofF-star expo- 
sures would result in non-zero residuals. From the 
amplitudes of the signals of the CO emission lines, 
measured at different directions from the star, and 
the noise level as shown in Figure 7, we estimate a 
signal-to-noise ratio of approximately 5—15. The 
noise arises partly from spurious mismatches in 
intensity between the on- and off-star spectra. 

To study the emission as a function of the angu- 
lar distance from the star, we divided our long-slit 
spectra into 79 sub-spectra (symmetrically around 
the maximum intensity representing the closest 
point to the star), one spectrum per pixel in the 
spatial direction. In this way we obtained 79 
spectra for every off-star slit position, each corre- 
sponding to 0.2" on the sky, if seeing is neglected. 
For every spectrum, the intensity of the three CO 
vibration-rotation emission lines [R(l), R(2) and 
R(3)] from the wind are measured. These spectra 
provide data representing a distance range from 
2" to a maximum of 7" away from the star, every 
slit position giving two series of data with three 
line-fluxes per spectrum, see Figures 8, 9, 10, and 
11; for example, the west position will sample the 
south-west and the north-west regions of the wind. 
The two sequences in every panel in the Figures 
represent these two series. In Figure 8 the decline 
of the added R(l), R(2) and R(3) intensities for 
the four slit positions are shown. 

Table 1 gives the observed intensities of the CO 
emission lines from the spectra representing the 
closest points to the star, i.e. (2 ± 0.5)" away 
from the star, as well as the mean intensity ratios, 
R(l)/R(3) and R(2)/R(3), for positions from 2.0" 
to 3.4". Since the S/N ratio at the observed spatial 
resolution is too low, the data will unfortunately 
not permit a useful plot of the intensity ratios as 
a function of distance; The scatter is too large. 
The best-fit slopes of the intensity as a function 
of angular distance from the star, dlog I/dlog/3, 
for the R(2) emission line are also given in the 
Table. The uncertainties quoted are pure mea- 
suring uncertainties, i.e. they do not include pos- 
sible systematic uncertainties. We obtained the 
best signal for our west position. The intensity 
decreases by a factor of about 40 from 2" to 7". 
A mean slope of the west, east, and north posi- 



tions is dlog//dlog/3 = —2.8 ±0.3. Here we omit 
the south measurement in order to lower the un- 
certainty. Naturally, the signal-to-noise ratios de- 
crease rapidly outwards. The uncertainty quoted 
is the measuring uncertainty. 

The R(l) emission line is situated close to the 
left (red) edge of the detector. The spectra show 
a peculiar rise in intensity here which makes the 
measured values of the R(l) lines subject to a sys- 
tematic uncertainty. Especially the east measure- 
ment of the R(l) line seems to be stronger than ex- 
pected when compared with the measurements in 
the other directions. The R(l) /R(3) ratios are also 
affected by this effect. The east and north spec- 
tra in Figure 6 show a spurious decline at lower 
wavenumbers (at the red end). This is due to a 
mismatch in intensity between the on- and off-star 
spectra at the edge of the detector array probably 
not caused by flat-fielding or other reduction pro- 
cedures. This spurious effect may, unfortunately, 
affect all frames. 

6. DISCUSSION 

We find a power-law dependence of the inten- 
sity as a function of angular distance on the sky, 
I oc Furthermore, in view of various uncer- 

tainties, e.g. in the positioning of the slit, we find 
that the measured emission line intensities in dif- 
ferent directions from the star are consistent with 
a symmetric wind. The wind is at least symmetric 
to within a factor of two in density. A true west 
and east position of 1.8" W and 2.2" E from the 
star would yield our measured values for a sym- 
metric wind. Corresponding values for the north 
and south positions would be 1.7" N and 2.3" S. 
Thus, a small error in the position of 0.2 — 0.3" 
(which is approximately the accuracy of the po- 
sitioning of the slit) in the measurement from a 
symmetric wind would yield the 'asymmetric' val- 
ues we measured. 

The relative intensities measured on the two 
nights can be checked at the cross-over points, cf. 
for example Figure 8. The intensities differ by ap- 
proximately 50%, probably reflecting the uncer- 
tainties in the slit positions, but also the accu- 
racy of the absolute flux calibration between the 
nights. The flux calibrations of the observations 
on the two nights are based on two different on- 
star measurements, which are subject to different 
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uncertainties, as discussed in Section 4. 

However, a comparison of the two direction in 
one slit position may reveal asymmetries in the 
wind. In Figures 9, 10, and 11, which represent 
the three lines measured, the increase in inten- 
sity at 3-4" in the 'north-east' data in the 'north'- 
panel is seen in all transitions. This is also the 
case for the 'south-east' data in the 'south'-panel 
which all show approximately the same morphol- 
ogy. Thus, this could indicate that there is an 
east-west asymmetry in the intensity measured, 
reflecting an asymmetric distribution of CO. This 
could be due to the additional bipolar outflow de- 
tected at larger scales in the aperture synthesis 
CO maps by Planesas et al. (1990b). Note that 
their highest resolution is approximately 6". 

Wc now make a simple analytic model of our 
observations by assuming a spherically symmetric 
and homogeneous wind with a constant mass-loss 
rate and a constant expansion velocity. We also as- 
sume the wind to be optically thin in the CO lines 
along the line-of-sight for rays from regions beyond 
a certain distance from the star. The adequacy of 
this latter assumption will be investigated below. 
Based on Eq. (5) by Gustafsson et al. (1997) the 
ratio of the wavelength-integrated, line-scattered 
intensity, Ico,i (ergs~^ cm~^ arcseconds^^), and 
the line-scattering flux f^ (ergs~^ cm~^ cm~^) as 
seen by the scattering molecules averaged across 
the line width but measured at the distance d, is 
found to be 



Ico,i(/?) 
fx 



206265 



32 Trie mn 



fu^i M X 



Ni(CO)/N(CO)-eco /IV 

f^v^d \f3j '^'> 

where /„<_; is the absorption oscillator strength of 
the line, eco is the fractional abundance of CO 
molecules, i.e. [CO]/[H], and Ni(CO) denotes the 
number density of CO molecules in the lower state, 
i, of the transition. Furthermore, fi is the mean 
molecular weight, d is the distance to the star, v^. 
is the terminal expansion wind velocity and /? is 
the angular distance from the star on the sky in 
seconds of arc. A feature of this general expres- 
sion is the minus third power dependence of the 
scattered intensity as a function of the impact pa- 
rameter on the sky, /?, which agrees within the 



uncertainties with our observation (cf. Figure 8). 
Eq. (1) will now be applied to the data for the R(2) 
line for which our data have the highest quality. 

The oscillator strength fu^i of the R(2) vibration- 
rotation CO-line at 4.6 /xm is 6.1 x 10~^ (Kirby- 
Docken and Liu 1978 provide oscillator strengths, 
which for our transitions are consistent with 
Goorvitch & Chackerian (1994) and Hure & 
Roueff (1996)3. See also Table 2). Thus, from 
Eq. (1) we find that the mass-loss rate as deduced 
from the emitted intensity from the R(2) line of 
CO, and, assuming n « 1.2, is 



M 



4.20 X 10"^° • 



Ni(CO)/N(CO) • eco 



Ico,i(/3) 

fA 



(Moyr-i) 



(2) 



where Vc is given in km s~^, d in parsecs. /3 in 
seconds of arc, Ico.i in ergs~^ cm^^ arcscconds^^, 
and f^ in ergs^^ cm^^ ^m^^. The fractional abun- 
dance of ^^CO molecules relative to hydrogen, 
eco. is assumed to be constant throughout the 
envelope and it is assumed that most oxygen is 
locked-up as CO molecules. From the literature 
we find fco = [COJ/pa] = 5 x lO""' (Knapp et al. 
1998), which means that eco = 2.5 x lO"'^. We 
now need to estimate the fraction Ni(CO)/N(CO) 
of CO molecules that are excited to the {v" = 
0, J" = 2)-level. We assume the population to 
be controlled by radiation transitions between the 
v" = and v" — 1 states. This radiation is 
originally supplied by the stellar photosphere at 
roughly the stellar effective temperature, T^s, but 
diluted by a factor of (R*/r)^ at a distance r from 
the star. Neglecting all loss or addition of photons 
in the optically thick spectral lines we find a char- 
acteristic temperature of the radiation from 



B{X,Tr) ^ B{X,ns) X (R*/r)2 



(3) 



This estimate gives a radiation temperature 
at a distance of 3" of about 340 K, suggest- 
ing an Ni=2(C0)/N(C0)-value of approximately 
4%. This temperature is several times higher 
than the kinetic temperature of the gas (the level 
is super-thermally excited). Detailed numerical 



'Note the missing cube in their Eq. (3) 
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simulations of the radiative transfer in an enve- 
lope model of o Ceti by Ryde & Schoier (2000) 
verifies the assumption of radiationally controlled 
populations and indicates population fractions, 
Ni=2(C0)/N(C0), that range from approximately 
3% to 10% at distances from 2" to 7". Here, 
we have adopted a value Ni=2(C0)/N(C0) of 5%, 
noting that a numerical modelling of the circum- 
stellar envelope will be necessary for a more de- 
tailed and accurate discussion. On the assumption 
that the terminal wind velocity Ve = 3kms~^, the 
mass-loss rate deduced from the measured inten- 
sity at 2" away from the star can be written as 

M = 0.10-Ico(/3)/fA (Moyr-i). (4) 

The observations at /3 « 2" of the east 
and west positions give Ico(/3)/fA = 3.1 x 
10~^/im(")^^ and of the south and north posi- 
tions give Ico(/3)/fA = 2.8 x 10" (")"^ where 
fA is estimated from the on-star spectrum. The 
mean mass-loss rate found by this analytic discus- 
sion is then M = 3.2 x lO"'^ Mq yr'^ for the east 
and west positions, and M = 2.8 x 1O~^M0 yr~-^ 
for the north and south positions. The uncertainty 
in the derived mass-loss rate (which is indepen- 
dent of the flux calibration since we use a ratio of 
fluxes in the derivation), are due to uncertainties 
in the position of the slit, the level population, 
the measured scattered intensity, and the on-star 
intensity. Considering these sources of uncertain- 
ties we estimate a conservative uncertainty in the 
mass-loss rate of typically a factor of 4. 

Our derived value of the mass-loss rate is in 
good agreement with values in the literature. For 
example. Young (1995) arrive at (3.6 ± 1.2) x 
10"^ Mq yr"i from the CO(J = 3-2) line for 
a distance of 100 pc and a eco = 3 x lO"**. The 
largest uncertainty in the estimates of mass-loss 
rates for miras has imtil recently been the distance. 
The Hipparcos satellite (ESA 1997) has radically 
improved this situation. Young's mass-loss rate 
corresponds to (2.8 ± 0.9) x IQ-'^M© yr'^ with 
our values of d and eco- Other estimates of the 
mass-loss rate lie between (1 — 5) x 10"'^ Mq yr~^, 
sec e.g. Danchi ct al. (1994), Loup et al. (1993), 
Planesas et al. (1990a), and Knapp & Morris 
(1985). It should be mentioned again that the 
radio profiles are difficult to model due to their 
asymmetry and require an additional component 



apart from for the standard CSE usually used. 
There may be a bipolar nebula seen pole-on (cf. 
the discussion in Knapp & Morris 1985). 

An assumption made in Eq. (1) is that the cir- 
cumstellar envelope is tangentially optically thin, 
i.e. along the line-of-sight, at distances from the 
star at which the measurements are made, i.e. 2- 
7". This optical depth is easily found from 



Ttang(p) = 



e'y fu^iM eco Ni{CO) 
meC^mH SX 
^ dl 



AjjLVe N{CO) 



2 ' 



(5) 



where p is the 'impact parameter' which denotes 
the closest distance from the star at the location 
on the sky of the observations and I is the variable 
of integration (tangential distance along the line- 
of-sight with the closest point to the star as the 
zero point). The integration should be calculated 
over a region of the envelope (— L < Z < L) where 
the differential Doppler shifts are smaller than the 
width of the line. Approximating the line profile 
with a step-function, we find this distance to be 



2L: 



d\ c 
p- . 

A Ve 



(6) 



An estimate of the tangential optical depth of the 
wind based on the derived mass-loss rate and the 
estimated expansion velocity, is obtained by solv- 
ing the integral (Eq. 5). We find that rtang ~ 7.6 x 
10^^ /p, where p should be in centimeters. Inwards 
of approximately 4" we find that Ttang(p)^l, i-e. 
that the wind is expected to be optically thick 
along a line-of-sight. Although the assumption 
that the optical depth is smaller than 1 is not 
fulfilled in the inner part of the observed region 
(/3 PS 2 — 7"), the CO emission declines smoothly 
as all the way in, indicating that the assump- 
tion of the wind being optically thin could nev- 
ertheless be correct. The dependence on angular 
distance of the scattered intensity is not expected 
to show the same behaviour in an optically thick 
and an optically thin wind. For the Ki scattering 
in the wind of R Scl (Gustafsson et al. 1997) the 
same situation is found. The /3~^-dependence of 
the intensity could be preserved throughout the 
observed parts of the wind if the circumstellar gas 
were inhomogeneous, which probably is a reason- 
able assumption, (cf. the discussion in Gustafsson 
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et al. 1997). This fact could lead to a lowering 
of the optical depth if the structures ('clumps') 
of the wind are of a suitable characteristic size 
and have a low filling factor in the linc-of-sight. 
Note, that the measured decline of the intensity 
could also be achieved if the lines are optically 
thick and the mass-loss rate is not uniform but 
vary with time. However, this seems to require a 
fine-tuning of the mass-loss rate with time which 
does not seem realistic. A numerical simulation of 
the Mira wind (Ryde & Schoier 2000) will discuss 
the optical depths in more detail. 

In applying Eq. (2) it was assumed that the 
line radiation at distance r from the center of the 
star is represented satisfactorily by the measured 
on-star flux f^, corrected by the distance factor 
(d/r)^. In practise, the observed on-star spectrum 
will also contain strong contributions from stellar 
light scattered in the inner envelope in non-radial 
directions, provided that the gas is not devoid 
of scatterers there. The velocity divergence then 
brings the photon out of the core of the line profile 
at points along the new direction and the chances 
for it to be re-scattered are minimal. Therefore, a 
considerable fraction of the photons are lost 'side- 
ways' from the envelope due to scattering already 
in its inner regions, see also the discussion in Ryde 
& Schoier (2000). 

7. CONCLUSIONS 

The successful detection of scattered photo- 
spheric light in the circumstellar CO vibration- 
rotation lines of Mira has made it possible to ex- 
plore the structure of its intermediate circumstel- 
lar envelope (at distances from the star of approx- 
imately 100 - 1000 stellar radii). 

The variation of the emission with angular dis- 
tance from the star (2" < /3 < 7") is found to 
roughly follow a f3~^ behavior, (3 being the angu- 
lar distance on the sky. Our data indicate that 
the mass-loss rate, when averaged over about 100 
years, was constant over a period of about 1200 
years, assuming that the wind is optically thin 
along the line-of-sight for the measured angular 
distances. 

Our observations provide absolute fluxes scat- 
tered in the circumstellar CO R(l), R(2), and R(3) 
lines, and provide corresponding line ratios. Based 
on a simple analytic model, the mass-loss rate is 



estimated to be approximately 3 x 10^"^ M0 yr"-'^ 
which, in view of the uncertainties, is compati- 
ble with earlier estimates in the literature. A nu- 
merical modelling of the wind will be necessary 
for a more detailed and accurate discussion of the 
circumstellar envelope. Such a modelling will be 
presented in a forthcoming paper (Ryde & Schoier 
2000). 

We have also found that the envelope is approx- 
imately spherically-symmetric to within a factor 
of two in density. Note, however, the asymmetries 
found by, for example, Planesas et al. (1990b). We 
have found indications that a similar asymmetry 
is also present in the CO vibration-rotation emis- 
sion intensity at a distance of approximately 3" 
from the star in the east direction as compared 
with the west direction. It is interesting to note 
that the companion of Mira is located eastward 
from the star at a position angle of (108.3 ± 0.1)° 
(Karovskaet al. 1997). 
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for valuable comments on the manuscript. We 
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Fig. 1. — The slit positions for our observations of 
Mira's envelope. While the spectrograph slit was 
30" long, the slit indicated in the Figure is 14 " 
long to indicate the regions in which we detected 
resonance-scattered CO emission. The underly- 
ing intensity map is taken from Planesas et al. 
(1990b). It shows the CO(J = 2-1) line emis- 
sion obtained with the IRAM 30 m telescope, in 
this case for 47.7kms~^, which is close to the stel- 
lar systemic velocity. The difference between two 
contour lines is 2 Kkms~^. The IRAM beam is 13 
" in diameter. 
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Fig. 2. — A low-resolution ISO spectrum (Rss 100) 
is shown in full line. The dashed line indicates the 
transmission profile of the Johnson M-band used 
in the absolute flux calibration of the observations. 
The width of the tall box denotes the observed 
spectral range of the Phoenix spectrometer. 
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Fig. 3. — The observed on-star spcctnim of Mira 
is shown as a full line. The superimposed dashed 
line represents the synthetic spectrum generated 
from a OS-MARCS model atmosphere for the stellar 
parameters Teg = 2400 K, logg = — 1.0 [cgs], and 
solar metallicity and assuming a stellar radius of 
3.8 • 10^"^ cm and a distance of 128 pc. The model 
fluxes are multiplied by a factor of 1.4 in order 
to match the observed spectrum. The agreement 
is astonishingly good, suggesting that the abso- 
lute flux calibration is reliable. The dominating 
spectral features are due to photospheric water va- 
por and CO. The telluric water and CO lines arc 
marked (®). Also, the absorption lines due to the 
photospheric CO 1-0 R(l), R(2), and R(3) tran- 
sitions are indicated in the Figure. Circumstellar 
absorption and emission is notable in the cores of 
the CO lines, especially R(3). 
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Fig. 4. — The observed scattered spectrum from 
an oflF-star observation including the R(2) line is 
presented as a dotted line. The solid line indicates 
a scaled on-star spectrum. From the Figure we 
can see that the small emission peak in the photo- 
spheric CO R(2) absorption line is located at the 
same wavcnumber as the off-star scattered spec- 
trum. This is the result of some residual circum- 
stellar emission in the on-star spectrum. Circum- 
stcllar absorption along the linc-of-sight through 
the circumstellar shell is present to the right (blue 
shifted) from the emission. The on-star emission 
and absorption components should nearly cancel 
in the measured scattered intensities, which are 
the differences of our off- and on-star spectra. 
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Fig. 5. — The spectra of o Ccti (Mira) with the 
abscissa showing the wavenumber scale of the ob- 
server. The off-star spectrum is shown by a dot- 
ted line and the appropriately scaled on-star spec- 
trum is shown by a solid line. The three cir- 
cumstellar CO vibration-rotation emission lines, 
which arc marked R(l), R(2), and R(3), 'fill in' 
the corresponding on-star lines which are in ab- 
sorption. To the right (blue) of these lines the tel- 
luric CO vibration-rotation R(l), R(2), and R(3) 
lines are seen. These telluric CO lines and the 
two major telluric water vapor lines at 2152.5 and 
2156.5 cm~^ are denoted with ®-symbols. The 
stellar CO lines are shifted out of the telluric ones. 
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Fig. 6. — The resulting circumstellar CO 
vibration-rotation flux from the wind of o Ceti. 
The intensity is integrated over the full slit area. 
A scaled on-star spectrum is subtracted from the 
off-star spectrum. In the four panels the spectra 
from the observations with the slit lying tangen- 
tially 2" W, E, N, and S, respectively, of the star 
are shown. Note the different ordinate scales. 
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Fig. 7. — A set of emission spectra from the 
west position is shown to indicate the signal-to- 
noise ratio at different distances from the star. 
The spectra arc shifted verticaUy with respect 
to each other. The top spectrum is measured 
closest to the star, i.e. at a distance of 2". 
The subsequent spectra are measured at 2.2", 
2.8", 3.6", and 4.5" away from the star, respec- 
tively. The intensity scale is given in units of 
10~^*erg s~^ cm~^ arcsec"^ /um~^. 
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Fig. 8. - The decline of the intensity of the 
circumstellar emission as a function of angular 
distance from the star. The line is a pure 
power law (I cx l3~^). The intensities of the three 
vibration-rotation lines [R(l), R.(2), and R(3)] are 
added in order to increase the signal-to-noise ra- 
tio. The two different symbols represent the two 
different directions along the long-slit, as indicated 
in the different panels. 
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Fig. 9. — The decline of the intensity of the cir- 
cumstellar emission as a function of angular dis- 
tance from the star, as in Figure 8. Here the 
scattered light of the CO R(l) vibration-rotation 
line is shown. 
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Fig. 10. — The decline of the intensity of the cir- 
cumstellar emission as a function of angular dis- 
tance from the star, as in Figure 8. Here the 
scattered light of the CO R(2) vibration-rotation 
line is shown. 
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Fig. 11. — The decline of the intensity of the cir- 
cumstellar emission as a function of angular dis- 
tance from the star, as in Figure 8. Here the 
scattered light of the CO R(3) vibration-rotation 
line is shown. 
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Table 1 

Measured intensities of the CO emission lines, mean intensity ratios and best-fit slopes 





West 


East 


North 


South 


Intensities [erj 


g s ^ cm ^ arcseconds ^] : 






Ir(i) 


3.3 X 10-12 


2.0 X 10-12 


2.3 X 10-12 


0.9 X 10-12 


Ir(2) 


2.6 X 10-12 


0.9 X 10-12 


1.5 X 10-12 


0.8 X 10-12 


Ir(3) 


2.0 X 10-12 


0.8 X 10-12 


1.3 X 10-12 


0.8 X 10-12 


Line ratios: 










R(l)/R(3) 


1.7 ±0.2 


2.1 ±0.6 


1.7±0.3 


1.2 ±0.2 


R(2)/R(3) 


1.3 ±0.1 


1.2 ±0.2 


1.2 ±0.2 


1.1 ±0.1 


Slopes of R(2), dlog7R(2)/dlog/3: 








-3.3 (NW) 


-2.5 (NE) 


-3.0 (NW) 


-2.3 (SW) 




-2.8 (SW) 


-2.5 (SE) 


-2.7 (NE) 


-3.2 (SE) 
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Table 2 

Transition probabilities for some CO lines 



Line (J" - J') 


wavenumber [cm ^] 






R(l) (1-2) 


2150.856 


6.791 X 10-6 


12.57 


R(2) (2-3) 


2154.596 


6.127 X 10"^ 


13.55 


R(3) (3-4) 


2158.300 


5.850 X 10-'5 


U.U 



25 



